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ABSTRACT 

An algorithm is presented for the computation of the topol- 
ogy of a non-reduced space curve defined as the intersection 
of two implicit algebraic surfaces. It computes a Piecewise 
Linear Structure (PLS) isotopic to the original space curve. 

The algorithm is designed to provide the exact result for 
all inputs. It's a symbolic-numeric algorithm based on sub- 
resultant computation. Simple algebraic criteria are given 
to certify the output of the algorithm. 

The algorithm uses only one projection of the non-reduced 
space curve augmented with adjacency information around 
some "particular points" of the space curve. 

The algorithm is implemented with the Mathemagix Com- 
puter Algebra System (CAS) using the SYNAPS library as 
a backend. 

Categories and Subject Descriptors 

1.1.4 [Symbolic and Algebraic Manipulation]: Applica- 
tions; 1.3.5 [Computer Methodologies]: Computer Graph- 
ics, Computational Geometry and Object Modeling-Geometric 
Algorithms 

General Terms 

Algorithms 

Keywords 

Algebraic Curves, Subresultants Sequence, Generic Condi- 
tions, Topology Computation, Sturm-Habicht Sequence, Ex- 
act Geometric Computation 

Introduction 

The problem of computing the topological graph of algebraic 
curves plays an important role in many applications such as 
plotting [13] and sectioning in Computer Aided Geometric 
Design [15], [16]. A wide literature exists on the computa- 
tion of the topology of plane curves f[8]. [7], [TO ] . [TT ] . [T2], [6] 
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and [H]). The problem of computing the topology of space 
curves has been less investigated. In [JJ, Alcazar and Sendra 
give a symbolic-numeric algorithm for reduced space curves 
using subresultant and GCD computations of approximated 
polynomials. If their approach gives good practical results 
however it doesn't give a rigorous proof that a sufficient pre- 
cision is selected for all inputs in the computation of GCD 
of approximated polynomials. In [TJJ, Owen, Rockwood and 
Alyn give a numerical algorithm for reduced space curve 
using subdivision method. Their algorithm has a good com- 
plexity but the topology around the singularities of the space 
curve is not certified. We also mention the work in [7J, 
where two projections of a reduced space curve are used, 
and where the connection algorithm is valid under genericity 
conditions. 

To our knowledge, the general problem of computing the 
topology of non-reduced space curves is not investigated 
in the algorithmic point of view despite its significance in 
the problem of computing the topology of a real algebraic 
surface. 

We present a certified algorithm that computes the topol- 
ogy of non-reduced algebraic space curves. We compute 
the topology of a plane projection of the space curve and 
then we lift the computed topology on the space. The 
topology of the projected curve is computed using a clas- 
sical sweeping algorithm (see [10], [8]). For the computation 
of the topology of a plane algebraic curve, we present an 
efficient generic test that certifies the output of the algo- 
rithm in [S] . 

For space curves, we introduce the notion of pseudo-generic 
position. A space curve is said to be in pseudo-generic po- 
sition with respect to the (x, j/)-plane if and only if almost 
every point of its projection on the (x, y)-plane has only 
one geometric inverse-image. A simple algebraic criterion 
is given to certify the pseudo-genericity of the position of a 
space curve. From a theoretical point of view, the use of the 
notion of curve in pseudo-generic position gives us a 
rational parametrization of the space curve. The use 
of this rational parametrization allows us to lift the topol- 
ogy computed after projection without any supplementary 
effort. From a practical point of view, the use of the rational 
parametrization of the space curve makes the lifting faster, 
avoiding numerical problems. 

We need to distinguish two kinds of singularities on the 
projected curve. A certified algorithm is given to do so. 
Unlike previous approaches, our algorithm uses only one 
projection of the space curve and works for non-reduced 
space curves. We therefore avoid the cost of the second 



projection used by previous approaches. 
In the next section we describe the fundamental algebraic 
tools that we use in this paper. In Section 2, we present 
our contribution to certify the algorithm for computing the 
topology of a plane algebraic curve. Our algorithm itself is 
introduced in Section 3. We report on our implementation 
and experiments in section 4. 

1. SUBRESULTANTS 

Let Pi,P 2 £ Q[X,Y,Z] and C R := {(x,y,z) € 1R 3 | 
Pi(x,y,z) = P2(x,y,z) = 0} be the intersection of the van- 
ishing sets of Pi and P 2 . Our curve analysis needs to com- 
pute a plane projection of Cr. Subresultant sequences are a 
suitable tool to do it. For the reader's convenience, we recall 
their definition and relevant properties. For all the results 
of this section, we refer to [3J, for proofs. 

Let A be a integral domain. Let P = X/T=o a iX l an d 
Q — Yli=o biX 1 be two polynomials with coefficients in A. 
We shall always assume a p 0, b q ^ and p ^ q. 

Let P r (A) be the set of polynomials in A[X] of degree 
not exceeding r, always, with the basis (as an A-module) 
1, X, . . . , X r . If r < 0, we set P r (A) = by convention, and 
we will identify an element S = so + . . . + s r X r of P r (A) 
with the row vector (so, • ■ • , s r ). 

Let k be an integer such that ^ k ^ q, and let ^k- 

P„_ fc _i(A) x P p _ fc _i(A) -» P p+ ,_ fe _i(A) 

be the A-linear map defined by *fc(J7, V) = PU + QV, with 
M k (P, Q) the (p + q - k) x (p + q - k) matrix of * fe . As we 
write vectors as row vectors, we have 

/ Oq ... Op \ 



V b ... b g J 

That is Mo(P, Q) is the classical Sylvester matrix associated 
to P, Q. To be coherent with the degree of polynomials, we 
will attach index i — 1 to the i column of Mk(P, Q), so the 
indices of the columns go from to p + q — k — 1. 

Definition 1 For j ^ p+q — k~ 1 andO ^ k ^ q, let sr^. j be 
the determinant of the submatrix of Mu (P, Q) formed by the 
last p + q — Ik — 1 columns, the column of index j and all the 
(p + q — 2k) rows. The polynomial Srfe(P, Q) = sr^.o + . . . + 
sr^fc X k is the k th sub-GCD of P and Q, and its leading term 
s?k,k {sometimes noted sr^) is the k th subresultant of P and 
Q. So, it follows that Stq(P,Q) — sro is the usual resultant 
of P and Q. 

Remark 1 

1. For k < j • ^ p + q — k — 1, we have STkj = 0, because it 
is the determinant of a matrix with two equal columns. 

2. If q < p, we have Sr 9 = (6 9 ) P_9_1 Q and sr, = {b q ) p - q . 

The following proposition will justify the name of sub- 
GCD given to the polynomial Sr^. 



Proposition 1 Let d be the degree of the GCD of P and 

Q (d is defined because A is an integral domain, so we may 
compute the GCD over the quotient field of A). Let k be an 
integer such that k ^ d. 

1. The following assertions are equivalent: 

a) k < d; 

b) Sr fe = 0; 

c) sr fe = 0. 

2. srd 7^ and Sr^ is the GCD of P and Q. 

Theorem 1 Fundamental property of subresultants 

The first polynomial Sr^ associated to P and Q with sr^. 7^ 
is the greatest common divisor of P and Q. 

We will often call (Sri)i the subresultant sequence associ- 
ated to P and Q and (sr;j )ij the sequence of their subresul- 
tants coefficients. We will denote by lcoefx(/) the leading 
coefficient of the polynomial / with respect to the variable 
X. 

Theorem 2 Specialization property of subresultants 

Let Pi, P2 G A[F, Z] and (Sii(Y, Z) )i be their subresultant 
sequence with respect to Z. Then for any a £ A with: 
deg z (P(Y,Z)) = deg z (P(a, Z)); 
deg z (Q(Y,Z)) = deg z (Q(a, Z)) , 

(Sii(a,Z) )i is the subresultant sequence of the polynomials 
P(a,Z) and Q{a,Z). 

2. TOPOLOGY OF A PLANE ALGEBRAIC 
CURVE 

Let / € Y] be a square free polynomial and 

C(/):={(a,/3)GP 2 ,/(a,/3) = 0} (1) 

be the real algebraic curve associated to /. We want to 
compute the topology of C(f). 

For curves in generic position, computing its critical fibers 
and one regular fiber between two critical ones is sufficient 
to obtain the topology using a sweeping algorithm (see [8]). 
But for a good computational behaviour, it is essential to 
certify the genericity of the position of the curve. 

We propose an effective test allowing to certify the com- 
putation and connection, in a deterministic way. This is an 
important tool in order to address the case of space curves. 

Now, let us introduce the definitions of generic position, 
critical, singular and regular points. 

Definition 2 Let f € Q[X, y] be a square free polynomial 
and C(f) — {(a, f3) £ TR 2 : f(a,f3) = 0} be the curve defined 
by f. A point (a, j3) £ C(f) is called: 

• a x-critical point if dy f{ct, ft) ~ 0, 

• a singular point if dxf(ct, (5) = 9y/(a, (3) = 0, 

• a regular point if dxf(a, /3) 7^ or <9y/(a, j3) 7^ 0. 

With these definitions we can describe the generic condi- 
tions required for plane curves. 



Definition 3 Let f £ QfX, Y] be a square free polynomial 
and C(f) = {(a,P) £ K 2 : f(a,f3) = 0} be the curve de- 
fined by f. Let Af x (a) := #{/3 £ IR, suc/i i/iat (a,/3) is a 
x-critical point ofC(f) } . C(f) is in generic position for the 
x-direction, if: 

1. Va £ C,A/;(a) ^ 1, 

2. There is no asymptotic direction of C(f) parallel to the 
y-axis. 

This notion of genericity also appears in [12] and in a 
slightly more restrictive form in [BJ. Previous approaches 
succeed if genericity conditions are satisfied, but they do 
not guarantee to reject the curve if they are not; i.e, it 
does not decide genericity. So for some input curves the 
computed topology might not be exact. 

A change of coordinates such that lcoefy(/) £ (Q* is suf- 
ficient to place C(f) in a position such that any asymptotic 
direction is not parallel to the y-axis. It remains to find an 
efficient way to verify the first condition. This follows from 
the next propositions. We refer to [8], for proofs. 

Proposition 2 Let / £ Q [X, Y] be a square free polyno- 
mial with lcoefy(/) £ (Q*, Resy(/, dyf) be the resultant 
with respect to Y of the polynomials f , dyf and { ax, . . . , a;} 
be the set of the roots of Kesy (/, dyf) in C 

Then C(f) is in generic position if and only if 
Vi £ {1, . . . , I}, gcd (f(oti, Y), dyfipti, Y)) has at most one 
root. 

Let / £ (Q [X, Y] be a square free polynomial with 
lcoefy(/) £ Q* and d := deg r (/). We denote by Su(X,Y) 
the i th subresultant polynomial of / and dyf and srij(X) 
the coefficient of Y 3 in Sr;(A", Y"). We define inductively 
the following polynomials: 

* (X) = ^o-o(X) 
° { ' gcdtsro.opasr^X))' 

V* £ {l,...,d-l}, = gcd($ l _ 1 (X),sr l , I (X)) and 

Proposition 3 

1. MX) = UTi(X) and Vi, j £ {1, . . . , d - 1}, i # j => 

g cd(r 1 (x)rr J (x)) = i; 

2. Let k € {1, . . . ,d - 1}, a £ C. r fe (a) =0i/ anrf onZj/ if 
g cd(f(a,Y),dyf(a,Y)) = Sr k (a,Y); 

3. {{a, 13) £R 2 :f(a,l3)=&rf(a,0) = 0} = U^{(a,/9) £ 
R 2 :r fc (a) = Sr fe (a,/3)=0}. 

In the following theorem, we give an effective and efficient 
algebraic test to certify the genericity of the position of a 
curve with respect to a given direction. 

Theorem 3 Let f £ Q [X, Y] be a square free polynomial 

such that degy(f) — d, lcoefy(f) £ (Q*. Then C(f) is in 

generic position for the projection on the x axis if and only 

ifMk £ Vi £ {0,...,k-l}, 

k(k-i) ST k ,i(X)sT k k (X) - (i + l)sr fe , fc _i(X)sr fe:i+ i(X) = 

0modr fc (X). 



Proof. Assume that C(f) is in generic position and let 
a £ C be a root of Tk(X). According to Proposition [3] (2.) 
gcd(/(o, r), 9r/(a, F)) = Sr fe (a, Y) = E * =0 s ^.,(«) yJ - 
According to Proposition [51 Sik(a,Y) has an only root 

= - k^HiUv so Srfc («- y ) = s ^Ma)(Y - p) k . 

Binomial Newton formula gives 

Sr k (a,Y) =sr k , k (a)(Y-/3) k = sr fc , fc (a) ^to (^H*)* - '**- 
So by identification VA; £ {1, . . . , d - 1}, Vi £ {0, . . . , fc - 1} 
and Va £ C such that T k (a) = 0, 

k(k - i) sr fe:i (a) sr feifc (a) - (i + 1) sr fc , fe _i(a) sr fe , i+ i(a) = 0. 

It is to say that Vfc £ {1, . . . , d - 1}, Vi £ {0, . . . , k - 1}, 
k(k - i)sr fc ,i(A)sr fe k (X) - (i + 1) sr fclii _i(X) sr M+ i(X) = 
Omodr^X). 

Conversely, let a be a root of T k (X) such that 

k(k - i) sr k ,i(a) sr fe , fc (a) - (i + 1) sr feifc _i(a) sr fc ,i +1 (a) = 0. 

With the same argument used in the first part of this proof 
we obtain 

gcd(f(a,Y),d Y f(a,Y)) = Si k (a,Y) 

= Ej= Brw(")^ (2) 
= sr k , k (a)(Y - (3) k 

with 

fcsr feifc ( a ) 

Then we conclude that gcd(/(a, Y), dyf (a, Y)) has only 
one distinct root and, according to Proposition [5] C(f) is in 
generic position. □ 

Remark 2 Theorem [3] shows that it is possible to check 
with certainty if a plane algebraic curve is in generic posi- 
tion or not. If not, we can put it in generic position by a 
basis change. 

In fact, it is well known that there is only a finite number 
of bad changes of coordinates of the form X := X + XY, 
Y := Y , such that if C(f) is not in generic position then the 
transformed curve remains in a non-generic position. This 
number of bad cases is bounded by (^) , where c is the num- 
ber of distinct x-critical points of C(f) [S]. 

3. TOPOLOGY OF IMPLICIT THREE DI- 
MENSIONAL ALGEBRAIC CURVES 

3.1 Description of the problem 

Let Pi, P 2 £ Q[X, Y, Z] and 

Cr := {(x,y,z) £1R 3 : Pi{x, y, z) = P 2 (x, y, z) = 0} (4) 

be the intersection of the surfaces defined by Pi = and 
P2 = 0. We assume that gcd(Pi,P2) = 1 so that is a 
space curve. Our goal is to analyze the geometry of Cr in 
the following sense: We want to compute a piecewise linear 
graph of E 3 isotopic to the original space curve. 
Our method allows to use a new sweeping algorithm using 
only one projection of the space curve. 

To make the lifting possible using only one projection, a 
new definition of generic position for space curves and an 
algebraic characterization of it are given. We will also need 
to distinguish the "apparent singularities" and the "real sin- 
gularities". A certified algorithm is given to distinguish 



these two kinds of singularities. 

For the lifting phase, using the new notion of curve in pseudo- 
generic position, we give an algorithm that computes a 
rational parametrization of the space curve. The use of this 
rational parametrization allows us to lift the topology of the 
projected curve without any supplementary computation. 

3.2 Genericity conditions for space curves 

Let II Z : (x,y,z) G M 3 ^ (x,y) G R 2 . We still denote 
n z = U z \ Cm . Let V = IL(C R ) C K 2 be the curve obtained 
by projection of Cr. 

We assume that deg z (Pi) = deg(Pi) and deg z (P2) = 
deg(P2) (by a basis change, these conditions are always satis- 
fied). Let h(X, Y) be the squarefree part of Res z (Pi, P2) G 
(Q[X, Y]. With the above notation and assumptions we have 
the following "geometric" equality, H z (Cjh) = C(h). 

Definition 4 [Pseudo-generic position] 

Let C c := {(x,y,z) G C 3 |Pi(z, y, z) = P 2 (x,y,z) = 0}. 
The curve Cr is in pseudo-generic position with respect to 
the (x,y)-plane if and only if almost every point ofH z (C<c) 
has only one geometric inverse-image, i.e. generically, if 
(a, 0) G n z (Cc), then U z 1 (ct,P) consists in one point pos- 
sibly multiple. 

Let m be the minimum of deg z (Pi) and deg z (P2). 

The following theorems give us an effective way to test if a 

curve is in pseudo-generic position or not. 

Theorem 4 Let (Srj(X,Y, Z)) je ^ m} be the subresultant 
sequence and (sij(X, y)) j6 r m j be the principal subresul- 
tant coefficient sequence. Let (Ai(X,Y)) ie ^ 1 ^^ m y be the 
sequence of <Q[X,Y] defined by the following relations 

. A (X, Y) = 1; Q (X, Y) = h(X, Y); 

• For i G {1, m}, 

&i(X, Y) = gcd(e i _i(X, Y), S n(X, Y)), 

Ai(X,Y) = q.\ x ^) ■ 

Fori G {l,...,m}, letC{A t ) := {(x,y) G M 2 |Ai(ar, y) = } 
andC(h) := {(x,y) G IL 2 \h(x, y) = 0} then 

m 

1. h(X,Y)= UMX,Y), 



2. Knowing that h(X, Y) = Y\Ai(X,Y), so it is clear that 

i = l 

m 

C(h) = UC(Ai). 

i=i 

3. Assume that Cr is in pseudo-generic position with re- 
spect to the (x, y)-plane. Let i G {!,..., m} and (a, /3) 
G C 2 such that sn(a,/3) / and Ai(a,/3) = 0. Then 
Ai(X,Y)= %-.tx?) ) =^9i-i(a,j9) =0. Knowing that 
6i_i(X,y) =*gcd(e < _ 2 (X ) y),sr i _i(A- ) y)) ) so it exists 
di,d 2 G Q[X, Y] such that 

&i- 2 (X, Y) = di{X, Y)O i - 1 (X, Y) and 
su-i(X,Y) = d 2 (X,y)ei_i(X,y). In this way, 
Gi-i(a,/3) = ==► 6,-2(0,/?) = and sri_i(a )( 3) = 
0. By using the same arguments, 6^-2(0, /3) = => 
6>;_3(a,/3) = and sr;_2(a,/3) = 0. By repeating the 
same argument, we show sr;_i(a,/3) = .. . = sro(a,/3) = 
0. Because sr,(a, (3) 7^ 0, then the fundamental theorem 
of subresultant gives 

gcd((P 1 (a,f3,Z),P 2 (a,[3, Z)) = Sn(a,P,Z) = 
YTj =0 sr M - j ( Q ' /^) Z l ~ 3 ■ Knowing that Cr is in pseudo- 
generic position with respect to the (x, j/)-plane and 
Ai(a,P) = then the polynomial Srj(a, /3, Z) has only 



one distinct root which can be written 



l(<*,/3) 



de- 



sn,i(a,/3) (Z + 



pending on the relation between coefficients and roots of 

m 

a polynomial. So Sii(a, j3, Z) = sr ,M-j ( a ) P)^^ 1 = 

3=0 

Conversely, assume that Vi G {l,...,m}, V(x,y) G C 2 
such that sri(x,j/) 7^ and Ai(x,j/) = 0, we have 

m 

Sn{x,y,Z) = J2 BT^ i - j (x,y)z i - i = 
3=0 

i,i-i( x ,y) 



z + 



Let be an irreductible 



component of U z (C<a). Then there exists i G {1, . . . , m} 
such that O C C(Aj). Let (a,/3) be a point of O, such 
that Ai(a,/3) = and sri(a, /3) 7^ 0. Now if we define 
, we obtain that Srj(a, j3, 7) = 0, then 



i(<*,/3) 



(a, /3, 7) is the only point of Cc with (a, /3) as projection. 
So Cr is in pseudo-generic position with respect to the 
(a;, y)-plane. 



2. C(h) = \JC(Ai), 

i=i 

3. Cr is in pseudo-generic position with respect to the (x,y)- 
plane if and only i/Vi G {1, . . . ,m},V(x,y) G <D 2 such 
that sn(x,y) 7^ and Ai(x,y) = 0, we have 



Sn(x, y, Z) = sr i , i (a; ! y) (Z + 



Proof. 1. By definition, Vi G {1, . . . , m}, 

Ai{X, y) = & q^ 1 ( < x^y) S> ■ So by a trivial induction 



e (x,y) 
e m (x,y)' 



deg z (Pi) = deg(Pi) and deg z (P 2 ) = deg(P 2 ) imply 
sr m (X,Y) G Q* (see RemarlQ]). 

So e m (x, y) = g cd(e m _ipf, y), sr m (x, y) )= 1, then 

m 

n Ai(x, y) = e (x, y) = h(x, y). 

i=l 



□ 



The following proposition is a corollary of the third result 
of the previous theorem. If Cr is in pseudo-generic po- 
sition with respect to the (i,y)-plane, it gives a rational 
parametrization for the regular points of Cr. 



Proposition 4 Assume that Cr is in pseudo-generic posi- 
tion with respect to the (x,y)-plane and let (ct,[3, 7) G Cr 
such that sii(a,P) 7^ and Ai(a,f3) = 0. Then, 



sr^ 



(a, 13) 



isr it i(a,/3) 



(5) 



Remark 3 By construction, the parametrization given in 
Proposition [4] is valid when sr^a, /?) 7^ 0. If sri,i(a,/3) = 
then either Aj(a,/3) = for some j > i or (a,/3) is a a> 
critical point of C(Aj) (see section 3.3). 



The following theorem gives an algebraic test to certify the 
pseudo-genericity of the position of a space curve with re- 
spect to a given plane. 

Theorem 5 Let (Sr., (X, Y, Z)).^^ Q m j be the subresultants 
sequence associated to Pi(X, Y, Z) and Pz(X,Y, Z) and 
(Ai(X, y))iG{i,..., m } be the sequence o/ (Q[X, Y] previously 
defined. The curve Cr is in pseudo-generic position with 
respect to the (x,y) -plane if and only if 

Vie{l,... J m-i},Vje{0,...,i-i}, 

*(» - j) sr*,j {X, Y) sn,i(X, Y) - (j + 1) m,i-i{X, Y) 
sn,i+i{X,Y) = 0modA»(X,y). 

Proof. Assume Cr be in pseudo-generic position. Let 
i € {1, . . . , m - 1}, j € {0, . . . ,i - 1}, (a, (3) £ R 2 such that 
Ai(a,j9)=0. 

If 8Ti t i(a,P) = 0, then sri,i_i(a,/?) = 0, consequently i(j + 
1) srij+i(a, /3) sri,i(a, /?) — (i— j) sri,i_i(a,/3) sri,j(a, /?) = 0. 
If STi t i(a, j3) ^ 0, then according to Theorem [4] (3.) 

Sxi{a,p,Z) = £ Br li , w (a,/3)Z i -^ = 

sr» * (a, /3) f Z + sr --- l( "' / l ) , ) ' . Let 7 := - , then 

Sr«(a,/3,Z) = Esr,i,H(a^)^ H = sr M (a, /?) (Z - 7)' • 
3=0 

Using the binomial Newton formula we obtain Sri (a, /3, Z) — 

£ sr,i,i- j (a,j3)Z i - j = sr<,i(«,/9)E(j)(-7) W ^'- So by 
3=0 3=0 
identification, it comes that 

Vie{l,... 1 m-l},Vje{0,... l i-l} ) 

i(i-j') sr i]3 - (a, /3) ST iti (a, p)-(i+j) m,i-i(a, p) ST ilj+1 (a, (3) = 
0, V(a, /3), A; (a, 0) = 0. The reciprocal uses the same argu- 
ments. □ 

Remark 4 Theorem [5] shows that it is possible to check 
with certainty if a space algebraic curve is in pseudo-generic 
position or not. If it is not, we can put it in pseudo-generic 
position by a change of coordinates. In fact, there is only a 
finite number of bad changes of coordinates of the form 

X := X + AZ; Y := Y + ^Z; Z := Z, 
with A, /x £ (Q* such that if Cr is not in pseudo-generic po- 
sition then the transformed curve remains in a non-pseudo- 
generic position pp. 

Let us introduce the definitions of generic position, critical, 
singular, regular points, apparent singularity and real sin- 
gularity for a space algebraic curve. 

Definition 5 Let M(X,Y, Z) be the 2x3 Jacobian matrix 
with rows (<9xPi, <9i-Pi, c>zPi) and {dxP^.dyPi.dzPi). 

• A point p£ Cr is regular (or smooth) if the rank of 
M(p) is 2. 

• A point p€ Cr which is not regular is called singular. 

• A point p — (a,/3, 7) £ Cr is x-critical (or critical for 
the projection on the x-axis) if the curve Cr is tangent 
at this point to a plane parallel to the (y,z)-plane. The 
corresponding a is called a x-critical value. 




x 



Figure 1: Apparent and real singularities. 

Definition 6 [Apparent singularity, Real singularity] 

We call: 

1. Apparent singularities: the singularities of the projected 
curve V — I1 z (Cr) with at least two points as inverse- 
images (see figure]!]). 

2. Real singularities: the singularities of the projected curve 
T> — II^Cr) with exactly one point as inverse-image (see 
figure\V\j. 

Definition 7 [Generic position) 

The curve Cr is in generic position with respect to the 
(x,y)-plane if and only if 

1. Cr is in pseudo-generic position with respect to the (x, y)- 
plane, 

2. T> = IL(Cr) is in generic position (as a plane algebraic 
curve) with respect to the x-direction, 

3. any apparent singularity ofT> = II^Cr) is a node. 

This notion of genericity also appears in a slightly more re- 
strictive form in pp. 

The aim of the next section is to give an algorithm to 
certify the third point of the previous definition of generic 
position. We give also in this section an effective way to 
distinguish the real singularities from the apparent ones. 

3.3 Distinguish real singularities and appar- 
ent singularities 

In this section, we suppose that Cr is in pseudo-generic 
position and V — II z (Cr) is in generic position as a 

plane algebraic curve. 

Let (Tj(X))j e {i n } be the sequence of T polynomials as- 
sociated to the plane curve T> and (/J^X))^/^...^} be the 
sequence of associated rational parametrization (see ([3} )• 
Let (Srj(X, Y, Z))^ g r m j be the subresultant sequence as- 
sociated to Pi , P2 € <Q[X,Y,Z]. For any (k,i) £ {I, . . . ,m}x 
{0, . . . , k — 1} let, R kyi (X, Y) be the polynomial 



k(k-i) sr k>i {X, Y) sr fc k (X, Y)-(i+l) sr^-^X, Y) sr fc , l+1 (X, Y). 



Lemma 1 Let (a,b) £ K 2 such that sr kyk (a,b) / 0, the 
polynomial Srfc (a, b, Z) = XliLo 8T k,i(a, b)Z l £ IR[Z] has one 
and only one root if and only i/Vi £ {0, . . . , k— 1 } Rk,i(a, 6) = 
0. 

For any j € {1, . . . , n } we define the sequences 
(uk(X))ke{i,...,j} and (v k (X)) ke { 2 j) by 

Ui(X) :=gcd(r i (J0,sr 1 ,i(A:,/3 i (X))) I 

Ufc(X) := gcd(sr fc , fc (X,&(X)),Ufc_i(X)) 

t)fe(X) := quo(«fe_i(X),w fc (X)). 
For fee {2, . . . , j} and i 6 {0, fc — 1}, we define (wfc,i(A )) by 

W M (X) :=w*(X), 

ui fc><+ i(X) := gcd(i? fe:l (X,/3 J (X)),«; fc , l (X)). 

Theorem 6 For any j S {1, . . . , n }, let (Tj. k (X) ) fce{1 ... j} 
and (Xj,k(X) ) fee the sequences defined by the following re- 
lations 

T jA (X) = qno(T j (X),u 1 {X)) and T ilk {X) := w k , k (X). 

XjM x ) : = quo(TO fei0 (x),r 3 , fc (x)). 

1. For any root a of Tj }k (X), the x-critical fiber (a,ft(a)) 
contain only the point (a, ft (a), 7j(a)) t«£/i 7j(a) := 

~ ksT k , k ( a ,0 ■(<*)) ' so ( a >#n a )) ls a real singularity. 

2. For any root a of Xj.k(X), (a, ft (a)) is an apparent 
singularity. 

3. Cr is in generic position if and only if for any (j, k) £ 
{2,...,n}x{2,...j} Xj - fc (X) = l. " 

Proof. 1. Let a be a root of Yj. k (X) := w k>k (X) = 
gpd(R ktk -i(X,l3 j (X)),Wk,k-i(X))- Then w k ,k'-i(a) = 
R klk - 1 (a,l3j(a)) = 0. 

Wk,k-i(X) := gcd(R k . k ^2(X,f3j(X)),w ktk ^ 2 (X)), so 
w ktk -2(a) = R k , k -2(a,Pj(a)) = 0. 

By induction, using the same argument, it comes that for 
i from to (k — 1), w kt i(a) = R k ,i(a, ft(a)) = 0. 
Wk,o{X) := v k (X), so v k (a) — 0. Knowing that v k (X) := 
quo(ttfc_i(X), u k (X)); u k (X) and u k -\(X) are square 
free, then u k -i(a) = and u k (a) 7^ 0. Knowing that 
u k (X) = gcd(sr k , k (X,l3 3 (X)),u k - 1 (X)), then 
sr fcife (a,ft(a)) / 0. 

u k -i(X) = gcd(sT k - llk - 1 (X,/3j(X)),u k -2(X)) and 

Ufc_i(a) = 0, so sr fc _i l j,_i(a, y 9j(a)) = Ufc_ 2 (a) = 0. 

By induction, using the same argument, it comes that for 

i from to k — 1 sr^^a, ft (a)) = 0. 

For i from to — 1 sri,j(a, ft (a)) = and 

sr k , k (a, ft (a)) ^ 0, so by the fundamental theorem of 

subresultants, 

gcd(Pi(a, ft (a), Z), P 2 (a, ft (a), Z)) = Sr fc (a, ft (a), Z) 

= Y^i=o SI k,i( a ; Pi( a )) zl ■ Knowing that 

gcd(Pr(a, ft (a), Z), ft (a, ft (a), Z)) = Sr fc (a, ft (a), Z) 

= X^i=o srfe . i ( a 'A( Q: ))^ 1 an( ^ f° r * from to (k — 1), 
R k: i(a, ft (a?)) = then by the previous lemma the poly- 
nomial gcd(Pi (a, /3j (a) , Z) , Pi (a, ft (a) , Z) have only one 

root -vlV) — _!itiMWl 
root 7, to).- fexsrfc _ fc(aift{Q)) - 

2. Let a be a root of the polynomial 

XjAX) := quo(w M (X),r 3lfc (X)). Then w M (a) = 
and rj jfe (a) = w k , k (a) 7^ because Wfc,o(A) and Fj, fe (X) 
are square free. For i from to k — 1, knowing that 



w M+ i(X) := gcd(ft M (X,ft(X)),w M (X)), io fe ,o(a) = 
and Wfc,fe(a) / 0, then it exist i 6 {0, . . . , k — 1 } such that 
R k ,i(a, f3j (a)) ^ 0. So by the previous lemma the poly- 
nomial Sr fe (a, ft(a), Z) = £) i=0 sr fei j(a, ft(a))Z* has at 
least two distinct roots. 

By definition Wfe,o(-X") := wjfc(-X'), so v&(a) = 0. Knowing 
that Wfe(X) := quo(ti fe _i(A'),u fe (X)); u fe (X) and u k -i_(X) 
are squarefree, then Ufc_i(a) = and u k {a) / 0. 
Ufc-i(a) = 0, ?ife (a) / and 
Ufe(X) = gcd(srfe ife (X,ft(X)),w fc _i(X)) imply 
sr k , k (a, ft (a)) / 0. 

itfc_i(X) = gcd(srfe_i i fe_i(X,ft(X)),Ufc_ 2 (A)) and 

Ufe-i(a) = imply sr fc _i, fe _i(a, ft (a)) = u fc _ 2 (a) = 0. 

By induction, using the same argument it comes that for 

i from to (k — 1) sr 4ji (a, ft (a)) = 0. 

For i from to (k — 1) sri i i(a, ft (a)) = and 

srfe,fe(a, ft (q)) ^ 0, so by the fundamental theorem of 

subresultants 

gcd(Pi(a,ft(a),Z),ft(a,ft(a),Z)) = Sr fc (a,ft(a),Z) = 
Et srM(a,ft(a))Z\ 

gcd(P 1 (a,ft(a),Z),P 2 (a,ft(a),Z)) = Sr fc (a, ft (a), Z) 
and Srfc (a, ft (a), Z) has at least two distinct roots im- 
ply that (a, ft (a)) is an apparent singularity. 

3. Cr is in generic position if and only if any apparent sin- 
gularity of T> = IT z (Cir) is a node. Knowing that the 
apparent singularities of T> which are nodes are exactly 
those with a root of xi.2(X) as x-coordinate, so Cr is in 
generic position if and only if for any 
(j,k)€{2,...,n}x{2,...,j}, X j,k{X) = l. 



□ 

3.4 Lifting and connection phase 

In this section, we suppose that Cr is in generic posi- 
tion that means that Cr is in pseudo-generic position, 
V = IT z (Cr) is in generic position as a plane algebraic 
curve and any apparent singularity of T> = H z (Cn) is a 
node. 

To compute the topology of Cr we first compute the topol- 
ogy of its projection on the (x, j/)-plane and in second we lift 
the computed topology. 

As mentioned in section 2, to compute the topology of a 
plane algebraic curve in generic position, we need to com- 
pute its critical fibers and one regular fiber between two 
critical ones. So to obtain the topology of Cr we just need 
to lift the critical and regular fibers of T> = IT z (Cr). 

Here after we explain how this lifting can be done with- 
out any supplementary computation for the regular fibers 
and the real critical fibers. And for the special case of the 
apparent singular fibers, we present a new approach for the 
lifting and the connections. 

3.4.1 Lifting of the regular points ofV = II z (Cr) 
The lifting of the regular fibers of T> = IT z (Cr) is done by 

using the rational parametrizations given in Proposition [4] 

3.4.2 Lifting of the real singularities ofv = IT(Cr) 

The lifting of the real singularities of T> = IT (Cr) is done 
by using the rational parametrizations given by 1. of Theo- 
rem [6] 

3. 4. 3 Connection between real singularities and reg- 
ular points 



(a,b2,c 2 ) 



(a, ft 72) • 

I 



(a, (i) 




(a, ft 72) 



Figure 2: Connection between real singularities and 
regular points. 



For a space curve in pseudo-generic position, the connec- 
tions between real singularities and regular points are ex- 
actly those obtained on the projected curve using Grandine's 
sweeping algorithm [8] (see figure [2}. 

3.4.4 Lifting of the apparent singularities 

The lifting of the topology around an apparent singularity 
is a little more complex. Above an apparent singularity of 
T> — II z (Cr) we have firstly to compute the z-coordinates 
and secondly to decide which of the two branches pass over 
the other (see figure [3}. We solve these problems by analyz- 
ing the situation at an apparent singularity. 

According to Theoremg](2.), V = II z (Cr) = U C(A;), so 

i = l 

an apparent singularity is a cross point of a branch of C(A;) 
and a branch of C(Aj) with i,j G {1, . . . , m}. So we have 
the following proposition. 



Proposition 5 If 

such that Ai(a,/3) 



a, (3) is an apparent singularity of T> 
= Aj(a,P) = 0, then the degree of the 
polynomial gcd(Pi (a, P,Z), P 2 (a, P,Z)) G R[Z] will be (i + 
3). 

Let (a, /3) be an apparent singularity of V such that 
Ai(a,/3) — Aj(a,(3) — and 71, 72 the corresponding 2- 
coordinates. So by Proposition [5] and Proposition [1] 

sr ,o(a,/3) = ... = sr it i(a,/3) = ... = sr jtJ (a,/3) = = 

sr 1+ j_i, i+ j_i(a,/3) = 0. 

By Proposition [4] for any (a, b, c) G Cr such that Ai(a, b) 
and sr i,i(a,b) ^ we have c = — : 



Zi := 



-. So the func- 

z st i i (a,0J 

gives the z-coordinate of 



tion (x, y) . 

v ' * sr^ i{x,y) 

any (a, b, c) G Cr such that Ai(a, b) = and sr^^a, b) 7^ 0. 
Ai(a,/3) — but sri t i(a,/3) — 0, so the function Z% is not 
defined on (a,/3). The solution comes from the fact that 
the function Zi is continuously extensible on (a,/3). Let 
ui be the slope of the tangent line of C(A;) at (a,/3) and 
t G II* . Let 7i (t) := Z^ + tui) = - ^^"o - 
Knowing that the algebraic curve Cr hasn't any disconti- 
nuity, it comes lim t ^ + 7 ;(t) = lim t _ - ji(t) — 71. By 
the same arguments, if we denote U2 the slope of the tan- 
gent line of C(Aj) at (a,/3) and 7j(t) := Zj(a,f3 + tui) = 
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Figure 3: Lifting of an 


apparent singularity. 
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Figure 4: Connection above an apparent singularity. 



then lim t-*o+ 7j(<) = lim t ^o- 7j(t) = 72. 
The values ui, -1*2,71 and 72 are computed using Taylor for- 
mulas and certified numerical approximations. 

Now it remains to decide which of the two branches pass 
over the other. This problem is equivalent to the problem of 
deciding the connection around an apparent singularity. Let 
(a, 61, ci) and (0,62,02) the regular points that we have to 
connect to (a,/3,7i) and (a, ,0,72). The question is which 
of the points (a, 61, ci) and (0,62,02) will be connected to 
(a, [3, 71) and the other to (a, 72) (see figure [3])? In [I] 
Alcazar and Sendra give a solution using a second projection 
of the space curve but it costs a computation of a Sturm 
Habicht sequence of Pi and P2. Our solution does not use 
any supplementary computation. It comes from the fact that 
71 is associated to iti and 72 to U2- Knowing that tii is the 
slope of the tangent line of C(Ai) at (a, /3) and 112 the slope 
of the tangent line of C(Aj) at (a,/3), so (a,/3,7i) will be 
connected to (a, &i,ci) if (a, 61) is on the branch associated 
to u\. If (a, bi) is not on the branch associated to u%, then 
(a, 61) is on the branch associated to 112, so (a, /3, 72) will be 
connected to (a, &i,ci) (see figure HJ. 

Remark 5 For a curve in generic position any apparent 
singularity is a node, so the slopes at an apparent singularity 



Curve 


Pi(x,y,z) 


P2{x,y,z) 


Time (s) 


1 


2 — ; — 2 — ; — 2 ^ 

x + y + z — 1 


2 2 ; — 1 

x — y — z + 1 


0.032 


2 


2 — ; — 2 — ; — 2 

x + y + z — 1 


3 — ; — — 2 — ; — — 2 — ; — 1 — ; — 1 2 

or + 3x z + 3x2^ + z 6 + y* — xyz — yz^ 


0.659 


3 


(a; — 2y + 2z) + y + z — 1 


z — z — (x — 2y + 2z) + 3(x — 2y + 2z)y 


2.125 


4 


(x - 2y + 22)^ + + z l - 1 


y A - (x -2y + 2zf - {x - 2y + 2z)yz 


1.031 


5 


(x - y + zf + + - 1 


y z - (x - y + zY - (x-y + z)zY - z\(x - y + z) z + y*) 


1.6963 


6 


(a; - y + z)' z + y l + z 1 - 1 


((x~y + zY + y A + z'Y - 4((aj -y + zY + y A ) 


2.228 


7 


Or - y + zf + y 1 - 2(x - y + z) 


((x-y + zY + y' + z A Y - 4((x - y + zY + V z ) 


2.875 



Figure 5: Running time of experimentations. 
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Figure 6: Computed topology of curve 2 of table l5l 

are always distinct that is to say u\ 7^ 112- 

4. IMPLEMENTATION, EXPERIMENTS 

A preliminary implementation of our method has been 
written using the Computer Algebra System Mathemagix. 
Results are visualized using the AxeQ algebraic geometric 
modeler which allows the manipulation of geometric objects 
with algebraic representation such as implicit or parametric 
curves or surfaces. 

Since existing methods have no publicly available imple- 
mentations, table[5]only reports our experiments, performed 
on an Intel(R) Core machine clocked at 2GHz with 1GB 
RAM. 
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Figure 7: Computed topology of curve 7 of table [5l 



